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ABSTRACT 

We use the Sloan Digital Sky Survey to construct a sample of 625 brightest group and cluster 
galaxies (BCGs) together with control samples of non-BCGs matched in stellar mass, red- 
shift, and color. We investigate how the systematic properties of BCGs depend on stellar mass 
and on their privileged location near the cluster center The groups and clusters that we study 
are drawn from the C4 catalogue of Miller et al. (2005) but we have developed improved al- 
gorithms for identifying the BCG and for measuring the cluster velocity dispersion. Since the 
SDSS photometric pipeline tends to underestimate the luminosities of large galaxies in dense 
environments, we have developed a correction for this effect which can be readily applied 
to the pubUshed catalog data. We find that BCGs are larger and have higher velocity dis- 
persions than non-BCGs of the same stellar mass, which implies that BCGs contain a larger 
fraction of dark matter In contrast to non-BCGs, the dynamical mass-to-light ratio of BCGs 
does not vary as a function of galaxy luminosity. Hence BCGs lie on a different fundamen- 
tal plane than ordinary elliptical galaxies. BCGs also follow a steeper Faber-Jackson relation 
than non-BCGs, as suggested by models in which BCGs assemble via dissipationless mergers 
along preferentially radial orbits. We find tentative evidence that this steepening is stronger in 
more massive clusters. BCGs have similar mean stellar ages and metallicities to non-BCGs 
of the same mass, but they have somewhat higher a/Fe ratios, indicating that star formation 
may have occurred over a shorter timescale in the BCGs. Finally, we find that BCGs are more 
likely to host radio-loud active galactic nuclei than other galaxies of the same mass, but are 
less likely to host an optical AGN. The differences we find are more pronounced for the less 
massive BCGs, i.e. they are stronger at the galaxy group level. 

Key words: galaxies: clusters: general - galaxies: elliptical and lenticular, cD - galaxies: 
fundamental parameters 



1 INTRODUCTION 

The central galaxies in galaxy clusters seem to be special - in many 
cases, the differences are visually obvious, because central cluster 
galaxies often have extended envelopes (i.e. they are cD galaxies) 
and they are usually the brightest (and most massive) galaxies in 
their clusters. The term brightest cluster galaxy (BCG) has thus 
become synonymous with the term central galaxy. 

At first glance, it might seem evident that the location of the 
BCG at the bottom of the potential well of a cluster must be the 
cause for any property which distinguishes it from other (cluster) 
galaxies. However, BCGs are also the dominant population at the 
massive end of the galaxy luminosity function, and thus, their 
properties are influenced both by their large masses and by the 
cluster environment. It is very difficult to disentangle these two 
influences, because it is difficult to find equally massive non-BCGs 
for comparison. Since most BCGs are early-type galaxies, their 
properties are often compared with the known scaling relations for 
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elliptical galaxies. It has been claimed that BCGs lie on the same 
Fundamental Plane as other ellipticals jOegerle & Ho essel 199 1]), 
but that they lie off its projections (e.g. the Faber-Iackson and Ko- 
rmendy relations) in that they have lower ve locity dispersions and 
larger radii than predicted by these relations I Thuan & RomanishinI 



1981 



iHoessel et all 1 19871 : ISchomberj Il98l lOegerle & Hoessell 



199 11) . More recently, it has been claimed that the slopes of the 
Faber-Iackson and Kormendy relations chang e as a function of 
of galaxy luminosity for all elliptical galaxi es jLauer et alj 200i 
iDesroches et al.l l2006h . On the other hand, iBrough et alNlOO: 



find that the surface brightness profiles (and thus the radii) of 
BCGs depend on the host cluster properties. 

The formation mechanism of BCGs is also a subject of much 
debate. Early on, it was suggested that BCGs form when galax- 
ies sink to the bottom of the po tential well of a cluster and merge 
(term ed galactic cannibalism: lOstriker & Tremainelll975l : I White! 
1 19761) . However, iMerrittI ( Il985h argued that tidal stripping would 
reduce the masses of cluster galaxies to the point where dynamical 
friction is too slow for this to be a viable mechanism. These analy- 
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ses assumed that clusters are static entities. A further mechanism to 
form BCGs in situ in the cluster is star formation in cluster cooling 
flows. At the centers of clusters, gas reac hes high enough density t o 
cool and condense into the cluster core JsiiSl 19761 ; iFabiarjl 19941) . 
But while the mass deposition rates inferred from the X-ray lumi- 
nosities of cooling flow clusters are of the order of several hundreds 

I I ' 

to > lOOOMoyr" (e.g. lAllen e t al. 1996), observed st ar formation 

rates are at most ~ lOOMoyr"' (Crawford et al.ll 19991) . Recent X- 
ray studies have furthermore demonstrated that the clus ter gas does 
not cool below ~ 2 keV (e.g. IPeterson & Fabiar]|2006l) . Moreover, 
this scenario predicts that the stellar populations of BCGs should 
be young and blue, which is clearly inconsistent with observations. 

These scenarios were proposed before hierarchical struc- 
ture formation was fully established as the standard cosmological 
paradigm, and for simplicity they neglected many of the processes 
that ta ke place when clusters assemble through mergers. IPubinskil 
( 1 19981) used N-body simulations to show that a dominant galaxy 
forms naturally by merging of m assive galaxies when a clu ster 
collapses along filaments. Recentlv. IPe Lucia & Blaizo^ ( |2006|) in- 
vestigated the formation of BCGs in the context of the Millenium 
Run simulation ( Springel et al. 2005). In their model, the stars that 
make up BCGs today are formed in a number of galaxies at high 
redshifts, which subsequently merge to form larger systems. The 
final BCGs assemble rather late; by a redshift of z ~ 0.5, on aver- 
age about half of the final stellar mass lies in the largest progenitor 
galaxy. Since many of these mergers take place very late when most 
galaxies have converted the bulk of their gas into stars, the merging 
events are very nearly dissipationless and are not associated with 
significant star formation. This scenario is supported by observa- 
tions that demonstrate that BCGs exhibit little sc atter in luminos- 
ity over a wide range of redshifts liSandag&,1973; Schneider et al.l 



ll983l ; |Postman & Laueill995l ; lAragon-Salamanca et al.l 19981) , and 
that their color evolution is consistent with a passively evolving 
stell ar population that formed at h igh redshifts (zfo,m ~ 2 - 5). 

iBovlan-Kolchin et alj l l2006l) used two-component N-body 
simulations of galaxy mergers to show that the remnants of dissi- 
pationless mergers remain in the fundamental plane. They showed, 
however, that the locations of the remnants within the fundamental 
plane, and thus on projected relations such as the Faber-Jackson 
and size-luminosity relations, depend on the orbits of the merging 
galaxies. During cluster assembly, infall occurs primarily along 
filaments, suggesting that mer gers onto the BCG m a y take place 
preferentially on radial orbits. iBovlan-Kolchin et al] l l2006l) show 
that BCGs would then be predicted to lie on steeper Faber-Jackson 
and size-luminosity relations than field galaxies. 



Although BCGs are probably not formed in cooling flows, 
they are believed to play an important role in regulating the rate at 
which gas cools at the centers of groups and clusters. The central 
cluster galaxies often harbor radio-loud active galactic nuclei 
(AGN), which may provide the necessary heating to counteract 
radiative cooling. [Burns ( 1990) find that 10 out of 14 cD galaxies 
in cooling flow clusters are radio-loud, comp ared to 3 out of 13 
in clusters without cooling cores. However, iBest et ak I ( l2005bl) 
show that radio-loudness also depends strongly on other galaxy 
parameters such as stellar mass. 

This work aims at comparing the properties of BCGs to those 
of "normal" galaxies of similar stellar mass, thereby disentangling 
the influences of mass and environment. Our work is based on data 
from the Sloan Digital Sky Survey (SDSS; Sect.|2j. It has recently 
been noted that the SDSS photometric pipeline tends to underesti- 



mate the luminosities of BCGs - we present a method to retrieve 
more accurate magnitudes from the catalogued data in Sect. [2] and 
Appendix|A] To constr uct a sample of BC Gs, we refine the C4 sam- 
ple of galaxy clusters ( lMilleretalJl2005l) in the SDSS by improv- 
ing the selection of the central galaxy (Sect.O. We also improve 
the C4 cluster velocity dispersion measurement by developing an 
iterative algorithm to simultaneously determine the redshift, the ve- 
locity dispersion, and the virial radius of each cluster (Sect. [3}. A 
detailed comparison of the measured BCG positions, redshifts, and 
velocity dispersions to the original C4 values is given in Sect.|4] In 
Sect.|5]we investigate the occurence of radio-loud AGN in BCGs. 
The results demonstrate that our choice of BCG is superior to the 
original choice of BCG in the C4 catalogue. The radio properties of 
the BCGs are inv estigated in more detail in an accompanying paper 
teest et alj|2006l) . as are the implications of these results for mod- 
els of cluster cooling flows. We compare BCGs to control samples 
of non-BCGs that are closely matched in stellar mass; we analyze 
their structural properties, their positions with respect to the funda- 
mental plane and its projections, their stellar populations, and their 
line emission in Sect.[6l A summary of our paper is given in Sect.[7] 

Our refined BCG/cluster sample, together with the wealth of 
information available from the SDSS, also provides a local compar- 
ison sample for optically selected high-r edshift clusters suc h as the 
ESO Distant Cluster Survey, EDisCS jWhite et al.|[2005h . While 
this work focuses on the BCG selection and on properties of the 
BCGs, another paper (von der Linden et al., in prep) will concen- 
trate on cluster properties such as substructure, mass estimates from 
velocity dispersion as compared to total light, and the properties of 
cluster galaxies as a function of their distance from the BCG. 

Unless otherwise noted, we assume a concordance cosmology 
with n„, = 0.3, Qa = 0.7 and Hq = 100/jkm/sAlpc, where h = 
0.7. 



2 DATA 



The S loan Digital Sky Survey jStoughton e7Zll2002l ; lYork et al.l 
I2OOOI) is a survey of about a quarter of the extragalactic sky, ob- 
taining photometry in five bands (iigriz) of more than 200 million 
objects and spectra of up to a million objects. The observations 
are carried out in drift-scan mode on a dedicated 2.5m telescope 
at Apache Point Observatory, with a large-array CCD camera that 
allows near-simultaneous photom etry. The i maging data is reduced 
by an automatic pipeline, photo jLupton et al. 2001), and various 
classes of objects are then classified for subsequent spectroscopy; 
those galaxies with 14.5 < m,- < 1 7.7 and 11 < 24.5 m ag arcsec"^ 
comprise the 'main galaxy sample ' ( IStraussetal]|2002l) . The spec- 
tra are obtained using a fiber-fed spectrograph on the same tele- 
scope. On each spectroscopic plate, which has a circular field of 
view of radius L49, 592 object fibers can be placed. Due to the fi- 
nite fiber size, any two fibers on the same plate need to be spaced at 
least 55" apart. The fiber allocation is performed by a tiling algo- 
rithm, whic h maximizes the number of objects that can be observed 
telanton et al. 2003). In the case of a "fiber collision" (i.e. two ob- 
jects that are closer than 55"), no preference of objects is given 
within the usual constraints. 

Our analysis is based on the fourth data release (DR4) of 
SDSS, whose main galaxy sample provides spectra for more than 
500000 galaxies. 
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2.1 Spectral analysis 

A multitude of physical properties have been derived for galaxies 
in the spectroscopic database via stellar population synthesis fitting 
and are publicly availabl^l The stellar continuum of each galaxy 
is modelled as a sum of template spectra generated from pop- 
ulation synthesis models jTremonti et alj |2004| ; iKauffmann et alj 
l2003ah . These fits also lead to measures of the stellar mass- 
to-light ratio, star formation histories, and mean stellar ages 
jKaufFmann et al. 2003a. b). After subtracting the stellar contin- 
uum, emission line fluxes can be a ccurately measured, allo wing 
studies of the star formation rates i Brinchmann et al. Il2004h and 
AGN activity ( iHeckman et alj2004lKaufl^mann et al.ll2003d) . 



2.2 Photometry and Stellar Masses 

It has recently been noted that the SDSS photometry systematically 
underestimates the lum inosities of nearby BCGs jBernardi et al.l 
l200d iLauer et alj [200d) . The problem arises because the level 
of sky background is overestimated both for large objects and in 
crowded fields. This is not only a problem for BCGs, but for all 
large galaxies, and the problem is worse in dense cluster environ- 
ments. Since any estimate of the stellar mass of a galaxy is derived 
from its luminosity, it is crucial for our analysis to avoid biases in 
the luminosity measurement. 

Apart from the local sky background measurement, which 
is estimated and applied for each galaxy, photo provides a global 
sky measurement estimated over a whole field. In AppendixlA] we 
argue that by adding up to 70% of the difference between local 
and global sky to the radial surface brightness profiles (as pro- 
vided by photo), more accurate photometry can be achieved. We 
test this p rocedure using ape r ture ph otometry of 35 BCGs from the 
survey of IPostman & LaueJ ( 1 19951) that are contained in the area 
of sky covered by the SDSS DR5. We also test our corrected mag- 
nitudes by comparing with the photometry provided by 2MASS 
jSkrutskie et al .1120061) . Our correction method (described in detail 
in Appendix [AJ has been applied to the BCGs as well as to about 
200,000 unique galaxies at z < O.I, which form the basis of our 
comparison samples for the BCGs. 

Since many BCGs are observed to have surface brightness 
profiles which do not follow a simple de Vaucouleurs profile 
jGonzalez et al .l2005l : lBernardi et al.l200l) . we choose not to use a 
magnitude measurement that assumes a certain profile, or is sensi- 
tive to the profile shape (this includes Petrosian magnitudes, which 
include about 80% of the light from galaxies with de Vaucouleurs 
brightness profiles, but almost 100% of the light from galax- 
ies with exponential light profiles). Instead, we measure isopho- 
tal magnitudes, defined as the light within the radius r-^^aii where 
the ID surface brightness profile reaches a surface brightness of 
(23-1- 10 log(I -l-z)) mag/n" in the r-band (the redshift term accounts 
for cosmological surface brightness dimming). This is a relatively 
bright isophote limit, chosen both to avoid residual uncertainties in 
the sky background subtraction (cf. Fig. lAlt as well as to exclude 
light from the cD envelope present in some BCGs (which is noti- 
cable typic ally at surface brightn ess levels one or two magnitudes 
fainter, cf. I Gonzalez et alj|200^ . We refer to these magnitudes as 
iso23 magnitudes. 

Other studies of BCGs have quantified their luminosities 
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within metric apertures (e.g. lPostman & Lauej|l99^ . These stud- 
ies were focused on cD galaxies, whereas our study includes BCGs 
of a much broader range in mass and size. Hence metric apertures 
would include very different fractions of light for BCGs at the ex- 
treme ends of the mass range, cf. Fig. [3] and any result based on 
them would be very difficult to interpret. It can be argued that the 
iso23 magnitudes are also dependent on the profile shape: For shal- 
lower brightness profiles, a larger fraction of the total light / mass 
is missed. But an estimate of the total light is unfortunately not 
feasible, due to residual uncertainties in the sky subtraction. We 
therefore adopt the iso23 magnitudes as the least biased and least 
problematic luminosity measurements. However, we have verified 
that our (qualitative) results do not change if we use Petrosian mag- 
nitudes as an attempt to measure total magnitudes. 

Stellar mass estimates for the BCGs and the comparison 
galaxies are derived from these luminosity me asurements using the 
KCORRECT algorithm telanton & Roweisll2006l) . which is also used 
to determine the k-corrections for our galaxies. BCGs that were 
not observed spectroscopically (see Sect. 13. U are assumed to have 
a redshift identical to the cluster redshift. Just as the luminosities 
should not be taken as an estimate of the total light, the quoted 
stellar masses are not an attempt to measure the total stellar mass. 
But since our analysis is based on comparing objects with similar 
masses and colors, this is not an issue. 



2.3 Radio catalog 

iBest et al. I( l2005ah identified the radio-emitting galaxies within the 
main spectroscopic sample of the SDSS DR2, by cross-comparing 
these galaxies with a combination of the National Radio Astron- 
omy Ob servatory (NRAO) V ery Large Array (VLA) Sky Survey 
(NVSS: ICoiidon et al.|[l998l) and the Faint I mages of the Radi o 
Sky at Twenty centimetres (FIRST) survey jBecker et al.lll99^ . 
They then used the optical properties of the galaxies to separate 
the radio-loud AGN from the radio-detected star-forming galax- 
ies. This work has now been extended to include the DR3 and DR4 
data (Best et al. in prep), and these results were used to identify 
those galaxies that are radio-loud AGN. 



3 SELECTION OF CLUSTERS AND BRIGHTEST 
CLUSTER GALAXIES 

The basis of ou r cluster sample is the C4 cluster catalog 
jMiller et alj2005l) . The C4 catalog is derived using the SDSS spec- 
troscopic sample and is currently available for Data Release 3. It 
identifies clusters in a parameter space of position, redshift, and 
color. The algorithm assumes that at least a fraction of the cluster 
galaxies form a color-magnitude-relation. It identifies galaxies in 
clustered regions, wit h neighbors of similar colors, as "C4 galax- 
ies" (see lMiller et al.| [2005i. for a detailed discussion of this selec- 
tion). From these C4 galaxies, it reconstructs the local density field, 
and identifies the C4 galaxies at the peaks of this density field as 
cluster centers (coined the mean galaxies). 

The C4 catalog identifies 1106 clusters within 0.02 < z < 

0. 16. In order to ensure that our clusters span a large angular ex- 
tent compared to the minimum distance between fibers (55 arcsec), 
we limit our cluster sample to z < 0. 1 . At this redshift, the magni- 
tude limit of the spectroscopic sample corresponds to ~ -20, 

1. e. slightly fainter than an L* galaxy. This cut results in a starting 
sample of 833 clusters. 
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3.1 Selection of the brightest cluster galaxy 

Our aim is to find the galaxy closest to the deepest point of the 
potential well of the cluster. In many rich clusters, this choice is 
obvious, and the central galaxy can easily be recognized as a cD 
elliptical galaxy by its extended envelope. Typically, this is also 
the Brightest Cluster Galaxy (BCG). However, in some clusters the 
central, dominant galaxy may not be the brightest galaxy. An ex- 
ample is the cluster C4_2003 H , shown in Fig. IBll The obvious 
central galaxy is SDSS J215729.42-074744.5 at the center of the 
image, but it is 0.3 mag fainter than SDSS J215701.71-075022.5, 
about 6'west-south-west of the former. We identify the former as 
the BCG (but concede that the term is a misleading nomenclature 
in this case). 

The C4 catalogue lists two galaxies for each cluster that could 
be considered the BCG: the mean galaxy (described above) and the 
brightest galaxy from the spectroscopic catalog within 500 /r' kpc 
of the position of the mean galaxy, four times the velocity disper- 
sion, and without strong Ha emission. However, due to the problem 
of fiber collisions, the true BCG is not included in the SDSS spec- 
troscopic data for about 30% of the clusters and is thus missed by 
the C4 algorithm. 

An earlier version of the C4 catalog tried to correct for this by 
selecting a brightest cluster galaxy based on the photometric cata- 
log. This object was selected to lie within 1 Mpc of the cluster 
center, and have a color compatible with the color-magnitude rela- 
tion of that cluster. Visual checks revealed, however, that this did 
not provide a reliable BCG: out of a subsamble of 128 clusters, 17 
of the photometric BCGs identified by C4 were stars misclassified 
as galaxies, and 36 were spiral galaxies (some of these located at 
the edge of the cluster). 

To identify the BCG for each cluster, we use the following 
procedure: 

(i) Based on the cluster redshift and velocity dispersion given by 
C4, we estimate the virial radius of the cluster: 



1.73- 



zh Mpc 



(1) 



1000 km s"' ^ai^ + a^il+zaf 

(see lFinn et alj[200^ . As C4 lists velocity dispersions within 
fixed physical radii (0.5, 1, 1.5, 2, 2.5 Mpc), we use the 
minimum non-zero value of these different values in this step, 
(ii) Within the projection of the larger of Rim and 0.5 Mpc around 
the mean galaxy, we select the two brightest galaxies that meet 
the following criteria: 

- The concentration index c = Rgo/R^o is larger than 2.5 
(where Rgo is the radius containing 90% of the petrosian 
flux measured in the /-band amd R^q is the radius containing 
50% of this flux), and f racDeV_r> 0.5 (this is a measure of 
the contribution of the de Vaucouleurs profile to the SDSS 
model r-magnitude). These cuts select galaxies likely to be 
early-types. 

- The color is compatible with that of the C4 mean galaxy to 
within A(u - g) < 0.6 (unless one of the u magnitudes has 
a large en'or estimate), A(g - r) < 0.5, A(r - /) < 0.4, and 



^ A note on the cluster IDs: the IDs used in the DR3 version of C4 are not 
identical to those in the DR2 version. Since the DR3 catalog was released 
only within the SDSS collaboration, the DR2 IDs are the "official' IDs, so 
whenever possible, we identify clusters by their DR2 ID (e.g. C4_2003). 
For those clusters without a DR2 ID. we use the DR3 ID, and denote these 
as e.g. C4.DR3.2004. 



A(r - < 0.4 (these are the dimensions of the color criteria 
originally used in the C4 algorithm to identify clustering in 
color space). 

- The flag TARGET.GALAXY has been set and the flag 
SATURATED is not set (these criteria allow us to identify stars 
that have been misclassified as galaxies, but they also apply 
to some low-redshift, bright galaxies). 

- If the object has spectroscopic data, we require that the red- 
shift is within Az < 0.01 from the cluster redshift. 

The brightest galaxy that meets these criteria is our initial BCG 
candidate. However, it is possible that these criteria select a 
foreground elliptical. Thus, if the second brightest galaxy is 
more than one magnitude fainter than the brightest, we also 
consider it as a BCG candidate. 

(iii) We then loosen some of these critieria: 

- we require only c > 2.3, 

- there is no constraint on A(m - g), 

- TARGET.GALAXY does not need to be set, and 
SATURATED can be set. 

Galaxies which meet this second set of criteria and are brighter 
than (both) the candidate(s) selected in the previous step enter 
the list of candidates. Unfortunately, misclassified stars enter 
our candidate lists at this stage. 

(iv) (a) If this procedure returns only one candidate, which agrees 

with the spectroscopic BCG given by C4, then this is au- 
tomatically considered the correct choice. This is the case 
for 242 clusters. 

(b) If this is not the case, then the BCG candidates (those given 
by C4 and those identified by our criteria) are inspected 
visually. For this purpose, we examine 2'.5 x 2'.5 color im- 
ages of the BCG candidates (provided by the DR4 Catalog 
Archive Server 

[CASB- 

These thumbnail images allow us 

to 

- identify cD galaxies by their extended envelope 

- identify stars misclassified as galaxies 

- identify obvious foreground ellipticals 

In the same step, a color image of the cluster (encompass- 
ing a field slightly larger than i?20o) can be viewed, with 
the BCG candidates marked. An impression can thus be 
gained of how the positions of the candidates relate to the 
C4 cluster center(s) (as given by the mean galaxy and the 
geometric cluster center, based on the mean of the posi- 
tions of the C4 cluster galaxies), visible galaxy overdensi- 
ties and other structures in the field. If there is more than 
one galaxy left in the candidate list, this color image al- 
lows us to choose the (brightest) one at the center of the 
galaxy overdensity identified by C4. In case such an over- 
density is not apparent, we choose the brightest elliptical 
galaxy in the vicinity of the the mean galaxy. 

For 472 clusters, the BCG can be identified fairly easily 
by means of these thumbnail images and the cluster image. 

(c) In case the previous step does not allow an unambiguous 
choice of the BCG, we enlarge the candidate list by adding 
galaxies within 2 Mpc that meet the criteria cited above. 
Along with this set of thumbnail images, we also inspect 
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a color image where galaxies within Az < 0.01 from the 
cluster redshift and 2 Mpc from the mean galaxy have been 
marked (see Fie. lBlb . It is this image that allows us to vi- 
sualize the clustering and to follow the galaxy distribution 
at the redshift in question. BCG candidates belonging to 
neighboring clusters can be identified. Conversely, if the 
C4 cluster corresponds to a weak number density fluctua- 
tion within another cluster, the BCG of this larger cluster 
is chosen. We identify the BCG of 54 clusters in this step, 
(d) The remaining 65 clusters require further scrutiny. For 
some of these clusters, this is necessary because the galaxy 
we identify as the central, dominant galaxy is not con- 
tained in the list of candidates. In other cases, it is evident 
from the cluster images that the cluster is in fact part of 
a larger cluster. In this step, we use the Finding Chart, 
Navigate, and Explore tools of the CAS website inter- 
actively. We start by marking galaxies within 2 Mpc and 
Az < 0.01 from the mean galaxy (as before), but then alter- 
ing the radius, center, and redshift range to gain an impres- 
sion of the clustering. We identify the center of the nearest 
galaxy overdensity and choose the most likely BCG within 
it. The color and brightness criteria that were previously 
used to determine the BCG candidate(s) are not applied in 
this step. 

At this stage, a given galaxy may have been identified as the BCG 
for more than one cluster. We then keep only the cluster whose 
mean galaxy is closest to the BCG. This rejects 101 clusters as 
being substructures of other clusters. 



3.2 Determination of the velocity dispersion and the virial 
radius 

The C4 algorithm measured velocity dispersions within fixed radii 
(see Sect. l4.3l for a more detailed discussion of the C4 velocity dis- 
persions). We prefer to measure the velocity dispersion o"v,ci within 
the virial radius R2Q0, which can be related to the velocity dispersion 
and cluster redshift Zcx via Eq. ([TJ. We thus developed an iterative 
algorithm to determine these three quantities: 

(i) From the catalog of galaxies described in Sect. 12. II we select 
those that lie within 2i?s,^rt from the BCG, where Rstmt is deter- 
mined using Eq. (TJ with a^ a given by the average value of the 
cluster velocity dispersions measured by C4. We also require 
the galaxies to lie within Az < 0.025 of the cluster redshift. 

(ii) For our first estimate of z^ii we take the redshift given by 
C4, z^j"*. For the first estimate of a^ ,.\ (and Rjoo), we use the 
median absolute deviation of the starting sample of galaxies 
with respect to . We also limit this first estimate to be less 
than 500 km/s, a step that is necessary to exclude surrounding 
large-scale structure. 

(iii) For each galaxy, we calculate its velocity within the cluster 
rest-frame: 



(iv) From those galaxies within ±3(Tv,ci of Zci, and within R200 from 
the BCG, we re-determine Zci, CTyd (and thus a lso Rim) using 
the biweight estimator from lBeers et al. ll fT990h . 

(v) Steps |(iii)l and |(iv)] are repeated until convergence is reached, 
i.e. subsequent iterations differ by less than 0.1% in Zci and 
CTv.ci. Galaxies are allowed to re-enter the sample. The galaxies 



contained in the final sample are considered the cluster galax- 
ies. 

It may happen that the iteration finds an oscillating solu- 
tion, i.e. iteration n + 2 yields the same solution (and the same 
galaxies) as iteration n. In this case, the algorithm modifies the 
input tTy ci to a random value between the two solutions and 
continues. Also, the convergence criterium is relaxed by factor 
of two every 200 iterations. A maximum of 1000 iterations is 
allowed. 

(vi) To estimate the error on the velocity dispersion, we draw 
10 000 bootstrap realizations from the cluster galaxies and cal- 
culate the velocity dispersion of each. We adopt 68% confi- 
dence intervals as our estimate of the ±lcr error . 

For 55 clusters, the algorithm does not converge or its sigma- and 
radius-clipping subsequently remove all (or all but one) of the 
galaxies from the starting sample. Since these systems cannot be 
considered bound, they are removed from the cluster list. 

At this point, all the clusters are inspected visually. We check 
a color image of the cluster with the cluster galaxies marked, as 
well as the redshift histogram. We mark clusters with the following 
indications that either the choice of BCG or the definition of cluster 
membership may be improved: 

• The redshift histogram does not justify the value of the veloc- 
ity dispersion. For clusters in very rich environments, the ve- 
locity dispersion can be overestimated if galaxies of neighbor- 
ing groups and filaments are included. Clusters with a velocity 
dispersion exceeding 1500 km/s are automatically flagged, but 
also those with two or more spikes in the redshift histogram. 
59 clusters are marked in this step. 

• Some of the cluster images show that the selected BCG is 
not at the center of the clustering. In others, there may be an- 
other bright elliptical present which could be the BCG, or the 
true BCG did not meet the criteria used to select BCG candi- 
dates (i.e. mostly late-type BCGs). 67 clusters are listed for 
re-investigation in this step. 

• Clusters in which the BCG redshift (if in the spectroscopic 
database) differs from the cluster redshift by more than Az = 
0.002 are also flagged (15 clusters). 

In this second round of visual checking, the choice of BCG or 
the sigma- and radius-clipping limits can be changed. The former 
is done if the previously selected BCG is clearly associated with a 
substructure of a larger cluster, but also if there is a possibly better 
BCG candidate that was previously missed. Note that many of these 
alternative BCGs are in fact foreground objects - whenever possi- 
ble, we retrieve redshift information on these objects from CAS 
or the Nasa/IPAC Extragalactic Database, NElfl For 21 clusters, 
the "new" BCG has not previously been selected as the BCG of 
another cluster. For 35 "clusters", the BCG has been attributed to 
another cluster, and so these systems are considered infall regions 
of these other clusters and discarded from the list. The sigma- and 
radius-clipping limits are changed for 53 clusters to avoid galaxies 
in nearby structures being included as cluster members (the sigma- 
clipping is typically changed from 3cr to 2.5cr or l.Qa - note that 
this affects only the choice of galaxies from which cry.d and zci are 
determined; cluster membership is still defined to be within 3(Tv,ci 
and i?20o)- For 67 clusters, nothing is changed. 

At this stage, we discard those systems that only contain 2 or 
3 galaxies within 3(Tv(.i and IRioo- This leaves 625 entries in our 

|http : //nedwww . ipac . caltech . edu/| 
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Figure 1. Histogram of the number of spectroscopic members within 
±3(Tv.ci and IS200 for our tinal cluster sample. Those clusters with > 100 
members are grouped into a single bin. The cluster with the most members 
is C4J3R3.3031 (Abell 2199) with 263 members. 



cluster sample. Fig. [T] shows a histogram of the number of spec- 
troscopic members for our cluster sample. The original C4 catalog 
considered only clusters with at least 8 galaxies within l/j^'Mpc 
and Az = 0.02. This is a much larger volume than we probe, and 
it is thus not surprising that we typically assign fewer galaxies to 
each cluster. 

Our final cluster sample spans a large range in velocity dis- 
persion, from groups of < 200 km/s to clusters of over 1000 km/s 
(Fig.|2l(. The positions of our BCGs as well as the cluster properties 
determined by our algorithm are listed in Table [T] Those systems 
with velocity dispersions < 300 - 400 km/s are more likely to be 
galaxy groups than galaxy clusters. We use the terms 'cluster' and 
'brightest cluster galaxy' loosely in this paper to refer to both clus- 
ters and groups. 

Fig-Elpresents a gallery of BCGs, sorted according to the ve- 
locity dispersion of the parent cluster (every 18th BCG is shown). 
The appearance of the BCG is certainly a function of (Xv ci, but it is 
not a monotonic one. While the BCGs of groups are mostly fairly 
isolated, rather spherical elliptical galaxies, the BCGs of more mas- 
sive systems are in general larger and more elongated, they often 
have a cD envelope and are surrounded by many satellite galaxies. 



4 COMPARISON TO C4 
4.1 Selected BCGs 

For 31 clusters in the final sample, the BCG we identified corre- 
sponds to both the C4 mean galaxy and the spectroscopic BCG 
identified by C4. 19 of these were classified automatically, as there 
was no other candidate in our list. 

In 343 clusters, our BCG is the same as the C4 spectroscopic 
BCG, but not the mean galaxy (183 automatically classified). 

In 41 cases, the BCG is the same as the C4 mean galaxy, but 



Table 1. Excerpt of our catalog of Brightest Cluster Galaxies for the C4 
cluster catalog, along with the cluster properties derived by our algorithm. 
The complete catalog is available electronically. Note that columns (1) and 
(2) refer to the C4 cluster ID (see footnote in Section[TTJ; column (7) gives 
the number of galaxies from which the cluster redshift and velocity disper- 
sion were determined. 



ID.2 


ID.3 


«BCG [°] 


6scG [°] 




o-y [km/s] 


Wgal 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1000 


1000 


202.5430 


-2.1050 


0.087 


648«lo 


35 


1001 


1001 


208.2767 


5.1497 


0.079 


746"„ 


82 


1002 


1002 


159.7776 


5.2098 


0.069 


m)%^ 


90 


1003 


1004 


184.4214 


3.6558 


0.077 


966'4o 


127 


1004 


1005 


149.7174 


1.0592 


0.081 




21 


1005 


1006 


191.3037 


1.8048 


0.048 


340"„ 


24 


1007 


1009 


177.4721 


5.7008 


0.075 


404«3 


27 


1009 


1011 


198.0566 


-0.9745 


0.085 


6311^5 


36 


1010 


1012 


192.0112 


-1.6528 


0.088 




14 


1011 


1013 


227.1073 


-0.2663 


0.091 




42 


- 


1014 


220.1785 


3.4654 


0.027 


45934, 


105 


1013 


1015 


203.0701 


1.2233 


0.079 


39759 
' -74 


10 


1014 


1016 


175.2992 


5.7348 


0.098 


660«, 


55 


1015 


1017 


182.5701 


5.3860 


0.077 


5965, 


41 


1016 


1018 


154.9344 


-0.6384 


0.093 


455I37 


17 


1017 


1019 


227.8480 


-0.0593 


0.091 


509'' 


36 


1018 


1020 


214.3980 


2.0532 


0.054 


6O55I53 


69 


1019 


1021 


195.7262 


3.3174 


0.071 


496'* 


24 


- 


1024 


199.8197 


-0.9954 


0.083 


579*' 


37 


1023 


1025 


153.4095 


-0.9254 


0.045 


790^25, 


66 


1341 


1026 


155.6325 


2.3608 


0.072 




26 


- 


1027 


191.9269 


-0.1373 


0.088 


1020*', [ 


55 


1027 


1028 


199.1357 


0.8702 


0.080 




16 


1030 


1030 


206.1357 


2.9541 


0.077 


■^^ -73 


29 


- 


1032 


218.4964 


3.7780 


0.029 


570^0 


76 


1032 


1033 


211.4731 


-1.2045 


0.054 


184«, 


4 


- 


1034 


165.7398 


7.6039 


0.072 




20 


1034 


1036 


192.3087 


-1.6874 


0.085 


771*3 
' ' -67 


64 


1036 


1038 


151.8861 


0.5942 


0.097 


550*' 


29 


- 


1039 


186.8781 


8.8246 


0.090 


846*3 

-66 


50 


1037 


1040 


213.6360 


1.7316 


0.054 


2994, 


16 


1038 


1041 


179.3707 


5.0891 


0.076 


678**, 


62 


1039 


1042 


228.8088 


4.3862 


0.098 


857**, 


53 


1040 


1043 


168.3339 


2.5467 


0.074 


403«, 


30 


1041 


1044 


194.6729 


-1.7615 


0.084 


"'-81 


60 


1044 


1047 


197.3295 


-1.6225 


0.083 


521*'* 


22 


1045 


1048 


205.5402 


2.2272 


0.077 


828^'gQ 


75 




1050 


206.1075 


2.1099 


0.072 


514*t, 


14 


1176 


1051 


189.7348 


6.1584 


0.074 


486^2 


19 


1048 


1053 


147.9551 


1.1339 


0.063 


346^* 


13 


1049 


1054 


188.7581 


1.7986 


0.079 


577*' 

-65 


35 


1051 


1057 


177.8878 


5.1015 


0.075 


251^2 


11 


1052 


1058 


195.7191 


-2.5164 


0.083 


74962 

-65 


68 




1059 


156.4666 


1.1906 


0.097 


478'" 


12 
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Figure 3. A gallery of BCGs. Each image is 200 kpc on the side. From left to right, top to bottom: 

C4J)R3.3351 (o-v.d = 159km/s), C4.2042 (tTv.d = 195km/s), C4J)R3.1343 (cr^.d = 210km/s), C4.3275 (cTv.d = 233km/s), C4.3087 (cTv.d = 256km/s), 
C4JDR3.3201 (cTv.d = 264km/s), C4.DR3.3106 (o-v.d = 283km/s), C4.1224 (cr^.d = 299km/s), C4.3206 (cTv.d = 312km/s), C4.2065 (cTv.d = 324km/s), 
C4JDR3.3272 (cTv.d = 329km/s), C4.DR3.1366 (o-y.d = 342km/s), C4J)R3.2140 (o-^.d = 355km/s), C4.3059 (cTv.d = 365km/s), C4J)R3.3386 (o-v.d = 
378km/s), C4J)R3.1355 (o-v.d = 392km/s), C4.3068 (o-v.d = 403km/s), C4.DR3.3034 (o-v.d = 410km/s), C4.1025 (o-v.d = 425km/s), C4.1226 (o-v,d = 
435km/s), C4J3R3.1360 (o-v,d = 448km/s), C4.3055 (o-v.d = 467km/s), C4.DR3.1356 (o-v.d = 484km/s), C4.1024 (o-v.d = 500km/s). C4.1076 (o-v,d = 
509km/s), C4.1191 (o-v.d = 519km/s), C4.1073 (o-v.d = 536km/s), C4.DR3.3105 (o-v.d = 556km/s), C4.3009 (o-v.d = 583km/s), C4J3R3.1275 (o-v,d = 
617km/s), C4J)R3.3027 (o-v.d = 670km/s), C4.1058 (o-v.d = 721km/s), C4.DR3.3084 (o-v.d = 781km/s), C4.DR3.3349 (o-v.d = 884km/s), 
C4.3002 (o-v,d = 1156 km/s). 
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Figure 2. Distribution of the redsliifts and velocity dispersions we measure 
for our final cluster sample. Error bars have been omitted for clarity. The 
marginal histograms show the distribution binned in redshift (lower panel) 
and in velocity dispersion (right panel). 



Figure 4. The differences in the positions of the C4 mean galaxy and our 
BCG, expressed in Mpc. Clusters in which the BCG is neither the mean 
galaxy nor the C4 spectroscopic BCG are shown as crosses, the other clus- 
ters are shown as filled circles. 



not the C4 spectroscopic BCG. These are predominantly small sys- 
tems, where the C4 spectroscopic BCG belongs to another system. 

The BCGs of 210 clusters correspond neither to the mean 
galaxy, nor to the C4 spectroscopic BCG. Of these, 141 (i.e. 23% 
of the whole sample) are not in the spectroscopic catalog. 

In Fig. |4] we compare the positions of our BCGs with those 
of the C4 mean galaxies. For the majority of the clusters, these 
two positions fall within 500 kpc of one another. 53 BCGs lie 
farther than 500 h^' kpc from the mean galaxy and of these, 21 
lie farther away than 1 Mpc. At these distances, the original cluster 
center is well in the outskirts of the structures we identify. 

An example is C4_DR3_1283 (see Fig.|B2ll, where the BCG 
and the original mean galaxy are separated by 2.8 Mpc, equivalent 
to 1.6i?200 according to the velocity dispersion we measure. The 
mean galaxy is a IO'^Mq, star-forming galaxy at z = 0.099, in a 
rich field that is likely to be an infall region of the cluster. The BCG 
we identified (a4x 10" Mq elongated elliptical with a cD envelope) 
is at the center of a cluster of 22 other galaxies. It is curious that this 
cluster was not picked up by the C4 algorithm. The comparatively 
high redshift of the cluster (z = 0.095) may possibly play a role in 
this. The cases of the other clusters with large separations between 
the BCG we identify and the C4 mean galaxy are similar, though 
less striking. 

4.2 Cluster redshift measurements 

The original C4 algorith m measures the c luster redshift using 
the biweight estimator of iBeers et al. I i l 19901) applied to all spec- 
troscopic members within an aperture of 1 /i"' Mpc from the 

^ 500 /i"' kpc is the radius within which the C4 algorithm identifies its 
spectroscopic BCG 



luminosity-weighted geometrical center of the cluster and within 
Az = 0.02 of the peak of the redshift histogram defined by these 
galaxies. While 1 /j"' Mpc is comparable to R20Q for a cluster with 
a velocity dispersion of ~ 600 km/s, the corresponding redshift in- 
terval from which we determine the redshift of such a cluster would 
be only Az = 0.006(1 -1- zcci)- 

In Fig.|5] we plot the relative differences between the two red- 
shift measurements as a function of the velocity dispersion cr"j 
measured by C4 within 1 Mpc. Our new redshift lies outside 
the Itr^^i redshift interval for only a few clusters , and only one lies 
outside the 3cr^*^ limit. The most notable outlier is C4JDR3_2163, 
with a velocity offset of about 2400 km/s. In our cluster sample, 
C4_DR3_2163 is a group of four galaxies at a redshift of Zd = 0.070 
and a velocity dispersion of 225^]Qikm/s- Its redshift histogram 
shows another spike of galaxies at a redshift of 0.082 which can be 
associated with C4_2124(see Fig. lB3l and Fig. IB4t . We thus con- 
clude that the C4 algorithm considered these two structures as a 
single cluster, whereas our algorithm was able to separate them. 

4.3 Velocity dispersion measurements 

The C4 catalog provides five measures of a cluster's velocity dis- 
persion, measured within 0.5, 1, 1.5, 2, and 2.5 /r' Mpc from the 
positional centroid measured on the sky. A first estimate for the ve- 
locity dispersion is made using the biweight estimator for all galax- 
ies within Az = 0.02 of the estimated cluster redshift. The final ve- 
locity dispersion (expressed in the observer's frame) is recomputed 
from galaxies within the redshift interval equal to ±4crJ;^|. 

Fig. [6] shows a comparison of our velocity dispersions to 
the C4 velocity dispersion within the radius that best corresponds 
to our estimate of the virial radius. At low velocity dispersions 
(o"v.ci 5 600km/s), the measurements agree well for many clusters. 
At higher velocity dispersions (as measured by C4) our algorithm 
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Figure 5. The differences between tlie C4 and our cluster redshift mea- 
surement (expressed as relative velocity) compared to the velocity disper- 
sion cr^^i measured by C4 within 1 /j"' Mpc (expressed in the cluster rest- 
frame). The dashed line indicates o"*^'', , and the dotted line 3(t^'*,. As in 

v,c] v,c] 

Fig. \4\ clusters in which the BCG is neither the mean galaxy nor the C4 
spectroscopic BCG are shown as crosses, and as filled circles otherwise. 



Figure 8. The radio-loud fraction of BCGs identified by C4 (green, filled 
circles), compared to those identified by our method (black, filled circles). 
We also show the radio-loud fractions derived only from clusters where the 
two choices differ: C4 BCGs are shown as blue, open circles, our BCGs are 
shown as red, open circles. Arrows indicate that the result is an upper limit 
only. The points are slightly offset in mass for clarity. 
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Figure 6. Comparison of the cluster velocity dispersions measured by C4 
and by our algorithm. We adopt the C4 velocity dispersion measured within 
the radius that best coiTesponds to the virial radius measured from our al- 
gorithm. Those C4 values that are larger than 1200 km/s are plotted at this 
value. 



obtains lower values for the majority of the clusters. As was pre- 
viously shown for C4_DR3_2163, this is mainly caused by the fact 
that our iterative algorithm separates neighboring groups/clusters 
better than C4. 



5 RADIO-LOUD AGN ACTIVITY OF BCGS 

It is known that BCGs often host radio-loud AGNs (e.g. iBurnj 
1 19901) . It has previously not been investigated whether this is sim- 
ply a consequence of the strong dependence of radio-AGN ac tivity 
on galaxy stellar mass (^adio-loud ^ 

Ml^ l lBest etal.|[2005bl) ), or 
whether this is a special property of BCGs. With our large sample 
of BCGs, it is possible to disentangle the mass dependence and the 
influence of the cluster environment. Fig. |7] compares the fraction 
of galaxies that are radio-loud for the BCG sample with the results 
found for all SDSS galaxies that overlap the NVSS and FIRST sur- 
veys. BCGs of all luminosities / masses are more likely to be radio- 
loud than other galaxies of the same luminosity / stellar mass. This 
enhancement ranges from a factor of 10 at ma sses of 5 x lO'^Mp 
to less than a factor of two above 4 x 10"Mr.. lBest et al.l ( l2005bl) 
have argued that radio-AGN activity is fuelled from the hot gas en- 
velopes of galaxies. In this scenario, groups and clusters provide an 
additional hot gas reservoir, which boosts the radio-AGN activity 
of the central galaxies. This result, and its implications for the cool- 
ing flow model, a re investigated in more detail in the accompanying 
paper (IBest^t al i2006i) . 

In this paper, we use the enhanced radio-AGN fraction of 
BCGs as a diagnostic for the reliability of the BCG selection, 
i.e. are our BCGs indeed better tracers of the bottom of the clus- 
ters' potential wells than the original C4 BCGs? In Fig. [8] we re- 
peat the previous analysis for C4 BCGs. At the highest mass bins 
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Figure 7. The fraction of galaxies that are radio loud AGN, as a function of luminosity (left) and stellar mass (right). Results ai'e plotted for all galaxies at 
2 < 0.1 (open triangles), and for the BCGs (solid circles). Galaxies are considered radio-loud if their 1.4 GHz radio luminosity is greater than 10^' W/Hz, and 
they are not classified as star-forming. 
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Figure 9. The same as Fig.[8] butforthe C4 mean galaxies. 



(> 10" 'Mg), the C4 BCGs have a similar radio-loud fraction to 
our BCGs, but at lower masses, the radio-loud fraction is lower 
than in ours. We also investigate the radio-loud fraction in clusters 
where our method and the C4 algorithm select different BCGs (252 
clusters). Here, the difference becomes even clearer: our BCGs 
have a much higher radio-loud fraction than those identified by 
C4. 



In Fig.|9]we repeat this analysis on the C4 mean galaxies. A 
similar trend as for the C4 BCGs is seen. 

These results can easily be explained by the difference in se- 
lection algorithm: in a cluster where the BCG is the most massive 
and brightest galaxy and has spectroscopic information available, 
C4 will correctly identify it, and hence the agreement is good in the 
high mass bins. However, if C4 misses the 'real' BCG, for example 
due to fiber collisions, it classifies a less massive galaxy (typically 
not at the bottom of the potential well) as a BCG. Since these galax- 
ies are normal cluster galaxies, their radio-loud fraction is lower 
than that of BCGs of equal mass. Hence, the C4 algorithm results 
in an underestimate of the radio-loud fraction at low masses. 



6 OPTICAL PROPERTIES OF BCGS 

Our large sample of BCGs and the extensive SDSS database allow 
us to compare the structural properties of BCGs with those of non- 
BCGs in order to distinguish the roles of mass and environment 
in governing their properties. In order to perform the comparison, 
we construct a comparison sample of "control" galaxies for each 
BCG by finding its three closest neighbors in a space spanned by 
(the logarithm) of galaxy stellar mass, redshift, and g-r color. The 
"redshift axis" of this space is scaled by a factor of five, so that a 
difference of 0. 1 in log M-,, corresponds to a redshift difference of 
0.02 , and a difference of 0. 1 in g - r. The matching is performed 
in order of decreasing BCG mass, and galaxies are not allowed to 
enter the comparison sample more than once. 

By matching in redshift, redshift-dependent aperture effects 
are avoided. The matching in g - r ensures similar stellar popula- 
tions and mass-to-light ratios in the BCGs and their controls, i.e. 
effectively, early-type BCGs are matched to early-type galaxies. 
Without the g-r matching, there are more late-type galaxies in 
the control sample than the BCG sample. But since our method of 
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selecting BCGs is somewhat biased towards selecting early-type 
galaxies over late-types, we cannot unambiguously deduce that 
BCGs are more likely to be early-types. 

The pool of galaxies from which the control sample is drawn 
consists of all galaxies in the DR4 spectroscopic catalog that have 
not been identified as a BCG in our sample. Yet, at the very mas- 
sive end (logMj^/Me > 11.5) there are not enough non-BCGs to 
provide three control galaxies per BCG. On the other hand, since 
we draw comparison galaxies from the full DR4 database, whereas 
the C4 catalog is based on DR3, the control sample is "contami- 
nated" by BCGs that failed to enter our sample, particularly for the 
most massive galaxies.]^ Hence, for a large part of the analysis, we 
restrict the sample to BCGs with logM«/Mo < 11.3; this avoids 
the problem of finding three suitable non-BCGs for the compari- 
son sample, and also purifies the comparison sample since at very 
high masses, a significant fraction of the comparison galaxies may 
themselves be BCGs. With these criteria, we construct two com- 
parison samples, one for the full set of BCGs, and the other for the 
subset of BCGs with spectroscopic information. 

In order to study scaling relations over a larger range in mass, 
we construct two more comparison samples (one drawn from all 
BCGs and one of them for BCGs with SDSS spectroscopy) with 
only one matched galaxy. Restricting to one comparison galaxy per 
BCG minimizes the problem of lack of comparison galaxies at the 
high mass end. For this matching, we restrict our analysis to only 
early-type BCGs and comparison galaxies by requiring Mg - M,- > 
0.75 and fracDeV_r > 0.8. 

Our four comparison samples are summarized below: 

CS3p: A comparison sample of three matching galaxies for BCGs 
with log M-i,/Mq < 1 1.3. The galaxies are matched in mass, 
redshift, and g - r color. 

CS3s: Like CS3p, but for BCGs contained in the spectroscopic 
database. 

CS Ip: A comparison sample of one matching galaxy for each BCG 
(with no upper mass limit). The sample is matched in mass, 
redshift, and g - r color, and restricted to only early-type 
BCGs and comparison galaxies (Mg - M,. > 0.75 and 
fracDeV.r > 0.8). 

CSls: Like CSlp, but for BCGs contained in the spectroscopic 
database. 

The first two samples are used to compare the distributions 
of physical parameters for BCGs and non-BCGs. The latter two 
samples are used to analyze early-type galaxy scaling relations, 
and to probe them to the highest masses. 

In Fig. [To] and Fig. [TT| we present the distributions of a va- 
riety of photometric and spectroscopic parameters for the BCGs 
and the comparison samples CS3p and CS3s. For each parameter, 
we list the decimal logarithm of 1 minus the Kolmogorov-Smimov 
confidence level at which the null hypothesis that the distributions 
are drawn from identical parent populations is rejected (i.e. a 99% 
probability that the distributions are different will have a value of 
-2). 

By construction, the BCGs and the comparison sample are 

' We have also attempted to clean the matched sample by applying the 
algorithm described in Sect. 13.21 to these galaxies, and considering those 
which are the brightest in structures of more than four galaxies to be pos- 
sible BCGs. About one third of the matched galaxies are such BCG can- 
didates. Basing our analysis on the remaining BCGs and the respective 
matched galaxies does not qualitatively alter our results. 



identical in stellar mass, redshift, and color. Because both stellar 
mass and color are the same, the distributions of luminosity are 
also equivalent. 

Note that our stellar masses are calculated using the kcorrect 
algorithm applied to isophotal magnitudes that have been cor- 
rected for sky-subtraction errors. These masses are not the same 
as the stellar masses es timated using the methods described in 
iKaufltmann et al.l ( l2003al) . which we compare in Fig.[TT] panel (c). 
For the latter, the mass-to-light ratio determined from the contin- 
uum spectrum is applied to the SDSS petrosian magnitude. Since 
these magnitudes underestimate the luminosity, the galaxy mass is 
underestimated as well. We also find that non-BCGs show stronger 
gradients between their fiber colors and their iso23 aperture colors 
(Fig.[lO] cf. Sect. l6.3t . Color gradients imply that the mass-to-light 
ratio varies across the galaxy and this is not accounted for when 
deriving stellar masses from the spectra. 

6.1 Structural parameters 

6.1.1 Radii and surface brightness 

In agreement with previous studies, we find that BCGs are larger 
(Fig. [TOl panels [f] and [g]) and have lower surface brightnesses 
than non-BCGs (Fig. [TO] panels [h] and [i]). The difference is 
more prominent for the inner characteristic radius R<,o (defined 
as the radius containing half the galaxy's light measured within 
the r = 23 mag/n" isophote) than for the outer isophotal radius 
i?iso23, within which we measure the luminosity of the galaxy. This 
is also evident in the distributions of the concentration parameter 
c' = RisoiilRiQ'- for a given i?iso23. a BCG has a larger ^50 than a 
non-BCG. This indicates that the light profiles of BCGs are sys- 
tematically different to those of non-BCGs. To first order, these 
differences c an be explained by BC Gs having shallower light pro- 
files. Indeed. iGonzalez et al.l ( f2005h and lBemardi et all ( l2006h find 
comparatively large Sersic indices (and thus shallow profiles) when 
fitting BCG light profiles with Sersic profiles. 



6.1.2 Size-luminosity relation 

The sizes and luminosities of elliptical galaxies ha ve been shown 
to obe y the scaling R^q oc L", with a ^ 0.6 (e.g. iBemardi et al.l 
l2003bl). However , at the massiv e end, this relation steepens 
l Laueretal.ll2006l) . iBemardi etal] ( I2OO6I) argue that BCGs have 
larger radii an d that this steepen i ng is c aused by an increasing frac- 
tion of BCGs. lOesroches et al.l ( I2OO6I) still find a steepening after 
removing the C4 BCGs from their sample of SDSS elliptical galax- 
ies and argue that the steepening is not solely attributable to 'con- 
tamination' from a population of galaxies with intriniscally larger 
radii (BCGs). 

The top right panel of Fig. [T2] demonstrates that BCGs are 
larger than non-BCGs at all luminosities or stellar masses. Sym- 
metric linear fits to the individual data points yield very similar 
exponents for the radius-luminosity relation for the BCG sample 
and the comparison sample: 



□ ^ r «.65±().Ui 

"SO.BCGs 1^ . 

n „ I 0.63±0.02 

nso.CSls L 

However, we also find that the relation displays curvature, i.e. it 
steepens with luminosity. This is shown in the lower right panel 
of Fig. [12] The range of expo nents we find {a ~ 0.5 - 0.7) is 
broadly consistent with that of iDesroches et all 1 I2OO6I) . We note 
that a is only slightly larger for the BCGs than the non-BCGs, 
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Figure 10. Distributions of a variety of photometric parameters for tlie sample of BCGs (red histograms) and the comparison sample, matched in stellar mass, 
redshift and g - r (black histograms). The ordinate of each plot shows the fraction of galaxies in a particular bin. In the top left corner of each panel we 
list the logarithm of the Kolmogorov-Smimov probability that the two distributions are drawn from an identical parent population. From left to right, top to 
bottom, the panels show: stellar mass; redshift; Mg - M, ; M, ; fracDeV_r; the inner characteristic radius S50, defined as the radius enclosing half the light 
measured within the isophotal radius Siso23; the r = 23mag/n" isophote radius R\so2i (within which our magnitudes are defined); concentration parameter 
c' = /?iso23/^5o; average surface brightness yUso within S50; the average surface brightness jUiso23 within Siso23; ^'^is ratio (from the flux- weighted second 
moments measured by photo); and the color gradient between the iso23 and the fiber apertures; . 



Table 2. The 16%, 50%, and 84% percentiles (left, middle, right columns respectively) of the distributions presented in Fig. llOl The values for the matched 
sample are listed in the top rows (black), those for the BCGs in the bottom rows (red). 



log(MWMo) z Mg-Mr M, 



10.84 


11.04 


11.19 


0.05 


0.07 


0.09 


0.83 


0.88 


0.94 


-22.62 


-22.24 


-21.76 


10.84 


11.04 


11.20 


0.05 


0.07 


0.09 


0.84 


0.88 


0.94 


-22.63 


-22.24 


-21.77 




fracDeV_r 






R50 / kpc 






Siso23 / kpc 






c' 




0.86 


1.00 


1.00 


3.14 


4.08 


5.41 


10.76 


13.79 


17.20 


2.92 


3.34 


3.69 


0.89 


1.00 


1.00 


3.24 


4.49 


5.89 


11.14 


14.47 


17.96 


2.84 


3.23 


3.65 
















hia 




(8- 


'•)i.so23 - (g - 


- '■)fiber 


18.31 


18.65 


19.09 


21.86 


22.04 


22.21 


0.76 


0.88 


0.95 


-0.07 


-0.03 


-0.00 


18.41 


18.82 


19.26 


21.93 


22.12 


22.32 


0.81 


0.90 


0.96 


-0.06 


-0.02 


0.00 
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Figure 11. As Fig. 1101 but for a comparison sample matched to the BCGs in the spectroscopic database; showing the distributions of various (mainly) 
spectroscopic parameters. The first panel demonstrates the match in stellar mass. The other panels show: stellar mass as extrapolated from the mass-to-light 
ratio derived from the continuum spectrum; velocity dispersion; strength of the 4000A-break; H(5a index; the metallicity index [Mg Fe]'; the alpha-to-iron 
index Mgb/<Fe); Ha line luminosity; HyS line luminosity; the Balmer decrement Ha/H^; the [Oin]5007 line luminosity; and the line ratio [Nn]/HQ' (a projection 
of the BPT diagram). 



Table 3. The 16%, 50%, and 84% percentiles (left, middle, right columns respectively) of the distributions presented in Fig. llll The values for the matched 
sample are listed in the top rows (black), those for the BCGs in the bottom rows (red). 



log(MWMo) 




log(M*/Mo)spcctra 


(T / km s 


A, (4000) 


10.82 


11.03 


11.18 


10.91 


11.17 


11.36 


187 239 279 


1.89 2.00 2.09 


10.82 


11.03 


11.19 


10.92 


11.17 


11.37 


203 246 288 


1.91 2.02 2.10 




H<5a 






[MgFe]' 




Mgb/<Fe> 


log(LHff/io) 


-2.57 


-1.88 


-0.80 


3.11 


3.35 


3.56 


1.38 1.62 1.85 


5.81 6.14 6.53 


-2.72 


-2.07 


-1.09 


3.10 


3.37 


3.58 


1.45 1.70 1.93 


5.75 6.07 6.47 




log(LH/j/Z,o) 










log(i-[Olll]/io) 




5.55 


5.80 


6.09 


1.96 


3.11 


4.32 


5.70 5.96 6.22 


-0.13 0.09 0.25 


5.53 


5.81 


6.07 


1.34 


2.67 


3.94 


5.69 5.92 6.17 


-0.19 0.07 0.26 
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Figure 12. Left panel: the size-luminosity relation for the BCGs (red) and for the compaiison sample CSlp (black) along with the best-fit linear relations 
(dashed lines). Top right panel: the median radii in bins of magnitude. The error bars represent the 68% confidence levels divided by the square root of the 
number of galaxies in the respective bin. Lower left panel: the variation of a, the exponent of the size-luminosity relation, as a function of magnitude A/^, 
determined from galaxies within M, ±1.0. 



even at the highest luminosi ties. We do no t find the significantly 
steeper rela t ions c laimed by ISernardi et al.l ( l200d . a = 0.92) and 
iLauer et al] jlOod . « = 1.18 ). Both of these samp les include both 
BCGs and non- BCGs - the iBemardi et all | |2006|) study uses the 
original C4 BCGs, which we have sh own to be contaminated by 
non-BCGs, and the iLauer et alj ( |2006|) study is based on galaxies 
with My < —21 and detectable core radii. Such contamination from 
non-BCGs is likely to be most important at lower luminosities, and 
will thus mimic a steeper slope. Even if we take this effect into ac- 
count, our data do not support very large values of a; if we fit a 

relation to the BCGs at M,. 23.5 and non-BCGs at M,. 23, 

we find a value of a of only 0.93. We speculate that one possible 
source for the discrepancy could be the different definitions of 
used in the different studies. 



6.1.3 Ellipticity 

We calculate the axial ratios of the galaxies in our sample (Fig. 1101 
panel [k]) from the flux- weighted second moments as measured by 
PHOTO, i.e. 



bja ■■ 



1 - yg + u- 

1 + ^IWTu^ 



(3) 



where Q and U are the Stokes parameters listed in the photo 
database. We choose this measurement since it is not based on 
fitting a particular model to the surface brightness profile of the 
galaxy and it is also not as sensitive to the sky subtraction as isopho- 
tal ellipticity measures. The majority of BCGs and non-BCGs are 
round, with axis ratios bja > 0.8. Both samples exhibit a tail to 
lower axial ratios, but this tail is more prominent for the non-BCGs. 
The median axial ratio is very similar for both samples (0.90 for 
the BCGs and 0.88 for the non-BCGs). This is qualtitatively con- 



sistent with t he results of iRvden et al. Ul993h . It should be noted, 
however, that lPorter et al.l ( Il99ll) find that the ellipticity of BCGs 
incr eases as a function o f the radius at which it is measured (see 
also lGonzalez et al.ll2005l) . 



6.2 Dynamical Structure 

6.2.1 Velocity dispersion 

As in previous studies of elli ptical galaxies that used spectra taken 
within fixed-sized apertures l |j0rgensen et al.|[l995l : iBemardi et al.l 
l2003ah . we correct the galaxy velocity dispersion to its expected 
value at one-eighth of the effective radius: 



''fiber 



(4) 



of the SDSS fiber (1.5"), and rso is the inner characteristic radius, 
measured in arcseconds. Strictly speaking, since the iso23 magni- 
tudes do not attempt to measure the total galaxy light, is not 
exactly the same as the effective radius, but because the correction 
does not scale very steeply with radius, this difference is negligi- 
ble. This correction also assumes a universal velocity dispersion 
pro file. While this seems to be applicable to most elliptical galax- 
ies l |j0rgensen et alj[T995l) , it has not yet been demonstrated that it 
also applies to BCGs. 

We find that BCGs have systematically larger velocity disper- 
sions than non-BCGs ((Fig. [TT] panel [c]; this also holds for the 
uncorrected velocity dispersions). 



© 2007 RAS, MNRAS Q0Q.mt24l 



How special are Brightest Group and Cluster Galaxies ? 15 




•^dUso / l^t.so log ("/ l<""/s) 



Figure 13. The ratio of dynamical mass to stellar mass (within R50) for 
the BCGs (red) and the comparison sample CS3s (black). We have assumed 
thatC2 = (1.65)-. 



Figure 14. Projection along the fundamental plane BCGs (red) and the 
galaxies of CSls (black). The dashed lines show the respective best tit for 
a, keeping h = 0.8. 



6.2.2 Dynamical mass 

The larger radii and higher velocity dispersions of BCGs imply that 
they have larger dynamical-to-stellar mass ratios than non-BCGs. 
The dynamical mass within i?5o can be derived via a projection of 
the scalar virial theorem onto observable quantities: 



M, 



dyn.50 



(5) 



where C2 depends on the profiles of both the dark matter and the 
luminous matter components. If the former follows an NFW pro- 
file (Navarro et al. 1997 ), and the latter aiHer nquist ( 1990,i) profile, 
then C2 = (1.65)2 IPadmanabhan et al.ll2004) . For calculating the 
dynamical mass, we assume that = (1.65)^ : 



Mdyn.5o = (1.65)-—— 



(6) 



We also assume th at the stellar mass within ^ 50 is 50% of the stellar 
mass within ^iso23 l lPadmanabhan et alj2004l ). We find that the ratio 
of dynamical mass to stellar mass is indeed considerably larger for 
BCGs (Fig. |13t . This difference is likely the consequence of the 
position of BCGs at or near the centers of galaxy clusters. As a 
result, there is a greater contribution from the dark matter halo to 
the dynamical mass of the BCG. 



6.2.3 The fundamental plane 

Early-type galaxies seem to be well described by a two-parameter 
set of equations, as is evidenced by the Fundamental Plane: they lie 
on a plane in a coordinate system defined by the logarithmic val- 
ues of velocity dispersion cry. effective radius, and average surfac e 
brightness within the effective radius jPiorgovski & Davislll987h . 
The plane is typically expressed as 



Re ^ (t"!-" 



(7) 



While there is agreement that b ^ 0.8, the parameter a is dependent 
on filter bands and may also be sensitive to a variety of selection 
effects and the precise definitions of cTy, Re and Ig. Typical values 
of a ~ 1.2-1.6 are quoted in the literat ure (see for example the 
compilation of observed FP coefficients in lBemardi et a l. 2003c). 

We assume that b = O.B. and plot \og{RioI°g^) as a function 
of logcTy in Fig. [14] for the early-type BCGs and the comparison 
sample CSls. In this diagram, the BCGs and the comparison galax- 
ies only diverge for galaxies with small radii and/or high surface 
brightness (i.e. these are not the cD galaxies). The difference is in 
the sense that the velocity dispersions of the BCGs are larger. A 
symmetric fit (fitexy from lPress et aljl992h yields 



"BCGs 

flcsis 



1.96 ±0.10, 
1.61 ±0.07 , 



For the comparison sample, the value of a lies close to the val- 
ues th at have been measured in the near-infr ared (e.g. Pahre et al] 
Il998l fl = 1.53±0.08), and also in the SPSS l lBernardi et alj2003"3. 
a = 1.49 ± 0.05). For the BCGs, a is significantly larger, indicating 
that BCGs do not lie on the same fundamental plane as "normal" 
ellipticals. It is interesting to note that it is predominantly the small, 
low velocity dispersion BCGs which deviate from the generic fun- 
damental plane. 

The fundamental plane relation is essentially an expression of 
the virial theorem. If we write 



Mj, = C[L and L = ItiI^qR^q . 
Eq. (|5j can be rewritten as 



R 



IttGc'i'm, 



(8) 



(9) 



dyn.SO 



The deviation of the observed fundamental plane from the theoret- 
ical one (a = 2 and b = 1) is referred to as the 'tilt' of the fun- 



damental plane. The tilt implies that ^ 



varies for different 
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Figure 15. The median luminosity to dynamical mass ratio of BCGs (red) 
and tlie comparison sample CSls (black) in bins of galaxy stellar mass. It is 
the variation of this ratio which determines the tilt of the fundamental plane. 



Figure 16. Top panel: the Faber- Jackson relation for BCGs (red) and for the 
comparison sample CSls. The dashed lines show the best linear fits to the 
relations. Bottom panel: the variation of fS with M^. i.e. ji determined from 
galaxies within M, ± 1.0 for BCGs meeting the early-type criteria (solid red 
line) and the corresponding CS 1 s sample (solid black line). The dashed line 
shows the subset of these BCGs located in clusters with d"v,c1 > 400km/s , 
the dash-dotted line the subset of BCGs in clusters with o"y_ci < 500km/s . 
The typical error bars are displayed on the left. 



elliptical galaxies. The proportionality constant ci is an expression 
of the stellar mass-to-light ratio and varies for different stellar pop- 
ulations. C2 depends on the profile shapes of both the luminous and 
the dark matter components. If ci were constant, elliptical galax- 
ies would be structurally homologous systems. There are contra- 
dictory results in the literature as to whether it is predominantly 
non-homology or variation in L/Mjy„ 50 that is responsible for the 
tilt of the fundamental plane. 

We are unable to distinguish non-homology from variation 
of Z./Mdyn,5o with our data. When calculating dynamical mass, we 
assume C2 = (1.65)^ (Eq. (6)), but we caution that this approach 
necessarily neglects effects from non-homology. 



In Fig. [TS] we investigate how the pre-factor — 



''1 Milyn.50 



T, T— varies as a function of stellar mass for BCGs compared to 

non-BCGs. The results show that the variation is much smaller for 
the BCGs. This is an affirmation of our previous result that BCGs 
lie on a different fundamental plane than non-BCGs. It also demon- 
strates that this result does not come from a few outliers, but applies 
to the majority of galaxies with M^, < 10" ^^ Mg. Since 



varies so little, the BCG fundamental plane is closer to the expec- 
tations of the virial theorem (a = 2 and b = \). 

Again, it is for low mass galaxies that BCGs differ most 
from non-BCGs. The similarity between BCGs and non-BCGs 
at high stellar masses implies that the process(es) which cause 
this ratio to be approximately constant for BCGs also apply to 
massive non-BCGs. Possibilities include assembly history (e.g. 
the i nfluence of the orbital eleme nts during dissipationless merg- 
ers, iBovlan-Kolchin et all |2006|) , and the fact that both BCGs 
and massive galaxies in general are found in dense environments 
jKauffmann et alj|2004l) . 



6.2.4 Faber-Jackson relation 

Several studies suggest that BCGs follow a different relation 
between luminosit y and velocity dispersion than less massive 
elliptical galaxies jOegerle & Hoesse 1 11991I : I Lauer et al ] I2OO6I: 
iBemardi et"al]|2006ir Parametrizing this relation as L oc cr'', the 
canonical value isjS = 4, as can be seen from Eq. l[8j and Eq. l|9]l, as- 
suming that — .l^* -r- is constant. Most measurements reported in 
the literature are consistent with f} =^ A (e.g. J jernard i et al.l2003br) . 
However, for sam ples of BCGs jOegerle & Hoesselll99lh and very 
massive galaxies JLauer et alj2006l) . it is found that/3 > 4 i.e. cr in- 
creases less steeply with luminosity than predicted by the standard 
Faber-Jackson relation. 

In Fig.[T6l we show the Faber-Jackson relation for the BCGs 
and CSls. Symmetric linear fits to each sample yield: 

, ^ ^5.32+0.37 

I ^ ^3.93+0.21 

We find a slope that is compatible with the standard L oc cr* rela- 
tion for non-BCGs and we confirm that cr rises less steeply with 
luminosity for BCGs.|2| 

In the bottom panel of Fig. [161 we investigate how fj changes 
with luminosity. We find that for BCGs, j8 is approximately con- 
stant (within the typical error bars) and has a value ~ 5.5. For non- 

^ We find a similar change in slope of the Faber-Jackson relation when 
using the /T-band luminosities of those BCGs with 2MASS photometry 
and a set of comparis on galaxies. This is contrary to recent claims by 
iBatcheldor et alj j2006l) . who use a much smaller sample of BCGs for their 
analysis. 
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BCGs,/3 varies from values ~ 3 at the low-luminosity end to values 
~ 4.5 at high luminosities (a similar range of albeit over a larger 
lumi nosity interval, as was found bv lPesroches et all2006l) . 

iBovlan-Kolchin et ai] ( l2006h find that in simulations of dissi- 
pationless mergers, j8 increases with the eccentricity of the merger 
orbit. They also argue that BCGs are expected to form through 
anisotropic merging due to the filamentary structure surrounding 
galaxy clusters. 

We test whether /? depends on cluster mass by splitting the 
BCG sample according to cluster velocity dispersion, and fitting f} 
separately for the two samples. We allow the sample to overlap in 
(Tv.ci to gain higher statistical significance. We obtain the following 
results: 

iBCGs(c^v,d < 500km/s) oc a-"2±o.46 
ific&s ( 0-v.ci > 400 km/s ) oc a- ^ 

The two values of p are just consistent with each other within 
the errors, and thus we cannot draw strong conclusions. Our re- 
sults indicate that p is larger for BCGs in more massive c lusters. 
If the scenario put forward bv lBovlan-Kolchin et al.l ( l2006h is cor- 
rect, this might imply that the merger orbit eccentricity increases 
with cluster mass. Another explanation might be that the number 
of (dissipationless) mergers is larger for BCGs in more massive 
clusters. 
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Figure 17. The median value of Mgb/{Fe) as a function of velocity disper- 
sion for the BCGs (red) and the comparison sample CSls (black). 



6.3 Stellar Populations 

The availability of measurements of spectral indices for galaxies in 
the spectroscopic SDSS catalog allows us to investigate the stel- 
lar populations in BCGs and non-BCGs. Both the distributions 
of the 4000A-bre£j <: (Fig. \n\ panel (d), measured as D„(4000), 
Balogh et al.' '1999*) and the H6a index (Fig. [TT] panel (e); 
Worthey & Ottaviani 1997) demonstrate that the stellar populations 
of the BCGs and the comparison galaxi es are old, as is generall y 
found for galaxies in this mass range jKauffmann et al.l l2003bl) . 
The metallicity, measur ed using the index [MgFe]' (Fig.[TTl panel 
(f); iThomas et al.ll2003l) is also typical for giant elliptical galaxies. 
There is a slight indication that the stellar populations of BCGs are 
slightly older (larger £)„(4000), lower HSa, and higher metallicity), 
but this is only significant for the HSa index. 

It should be noted that these measurements apply only to the 
galaxy light contained within the fiber, i.e. the inner 3", whereas 
the samples are matched in g - r color within the 23 mag/n" 
isophote. In panel (1) of Fig. [TO] we find that the color gradient be- 
tween the fiber aperture and the iso23 aperture is more prominent 
in the non-BCGs. This is a confirmation of previous results tha t 
color gradients in BCGs are weak or absent jGarim et al.lll997b . 
whi le non-BCG elliptical g alaxies are generally redder in the cen- 
ter lEaBaAeraetanioOSb. The presence of color gradients is typ- 
ically attributed to metallicity gradients jjames et al.ll2006l) , how- 
ever, we do not find evidence for different metallicities in BCGs. 

We use the ind ex Mgb/(Fe) as an indicator of the a/Fe ratio 
jThomas et ai]|2003l) and we find that BCGs have a systematically 
higher Mgb/(Fe) value than non-BCGs (Fig.[TT] panel [g]). How- 
ever, this index is known to correlate strongly with velocity disper- 
sion, so this result is not independent of our previous result that 
BCGs have systematically larger velocity dispersions. In Fig. [T7] 
we plot Mgb/(Fe) as a function of velocity dispersion for the BCGs 
and the comparison sample CSls. Except in the outermost bins, 
we do find systematically higher a/Fe ratios in the BCGs. Higher 
a/Fe ratios can be interpreted as an indication that star formation in 



the galaxy occurred over a shorter time-scale dCranato et al .120041) . 
The enhanced radio-AGN activity we find in BCGs (Sect.|5) may 
explain why star formation has been shut off on shorter timescales 
in the BCGs. 



6.4 Emission line properties 

We also investigate the strengths of the four emission lines Ha, 
Bj3, [Oiii]/I5007, and [Nn]/i6584 that are commonly used to clas- 
sify galaxies according to whether their em ission line luminosit y 
is driven by star formation or AGN activity ( [Baldwin et ai]|l98l[) . 
We limit our sample to galaxies with a signal-to-noise SNR > 3 
in the respective line measurement(s). For the individual line mea- 
surements Ha/ ¥1/3/ [Oni] / [Nii] , this holds for 56% / 42% / 64% 
/ 50% of the BCGs and 66% / 47% / 77% / 63% of the comparison 
sample (note that this applies to the CS3s sample). Requiring that 
SNR > 3 in all four bands simultaneously leaves only 30% of the 
BCGs and 40% of the comparison sample. These numbers already 
indicate that the emission lines in BCGs are in general weaker than 
in non-BCGs, a result which is further confirmed by the distribu- 
tions of Ho- and [Oiii] line luminosities shown in Fig. [TT] panels 
(h) and (k). We note that it is particularly the high-mass BCGs 
in which the emission line strength is suppressed compared to the 
comparison sample. 

Fig.|18|shows the BPT diagram of BCGs and the comparison 
sample for those galaxies which satisfy SNR > 3 in all four bands, 
i.e. 119 (out of 391) BCGs and 472 (out of 1173) non-BCGs. Of 
these BCGs, 7 (i.e. 6% of the line-emitting sample / 2% of the 
complete sample) are classified as star-forming, 83 (70% / 21%) 
as AGN, and 29 (24% / 7%) as composite. Of the non-BCGs, 29 
(6% / 2%) are star-forming, 364 (77% / 31%) AGN, and 79 (16% 
/ 6%) composite. Our sample is too small to draw detailed conclu- 
sions from these numbers, except that for both samples, the emis- 
sion line flux is dominated by AGN-like emission. Their low [Oiii] 
luminosities place the BCGs somewhat lower in the BPT diagram 
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Figure 18. BCGs (red) and the comparison sample CS3s (black) placed 
into the BPT diagram. Galaxies to the left of the blue line are classified as 
purely star-forming {Kauflmann et al. 2003c ), galaxies to the right of the 
green line as purely AGN tKewley et aL2001,) . and galaxies inbetween the 
lines as composite. Only galaxies with SNR > 3 in all four line measure- 
ments are shown. 
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Figure 19. The Balmer break D„(4000) vs. M„ - Mg for the BCGs (red) 
and the comparison sample CS3s (black). Typical eiTor bars are shown in 
the lower right corner. The BCGs which were identified as star-forming 
from the BPT diagram are marked by circles. Note that only one of these 
is an early-type galaxy and is also in the CS 1 s sample. In the latter sample, 
there is no BCG and only one comparison galaxy with D„(4000) < 1.6. 



than the non-BCGs, i.e. the BCGs are almost exclusively classified 
as LINERs, whereas a few galaxies in the comparison sample could 
be classified as low-luminosity Seyferts. 

For case B recombination, the unattenuated value of the 
Balmer decrement is ~ 3 (more specifically, it is Hty /HS = 2.86 
in sta r-forming galaxies and Ha/HyS = 3.1 in AGN, lOsterbrockl 
Il989h . Remarkably, a considerable fraction of the BCGs in our 
sample have Balmer de crements below this value (Fig. [TT| panel 
[j]). It has been noted bv lKewlev etall | |2006|) that 33% of LINERs 
in the SDSS sample have Ylafttf} < 2.86. A possible explanation 
is that the fits to the stellar continuum are not entirely reliable for 
the most massive galaxies, which tend to have very strong metallic 
absorption lines in their spectra. 

6.5 Star formation in BCGs 

The results of the previous section suggest that BCGs do not 
have increased amounts of star formation with respect to the com- 
parison sample. This is somewhat surprising, as there are sev- 
eral BCGs known with signs of recent star f ormation (blue col- 
ors, Hff emission; see ICrawfordet"ai]|l999h , and some cluster 
cores are known to have extended Hce emission structures (e.g. 
ICrawford et al.ll2005l) . The latter are found t o occur exclusively in 
cooling core clusters ( ICrawfordet"ai]|l999l) . However, since the 
SDSS fiber only probes the inner galaxy, it is inse nsitive to such 
surrou nding filaments. As for nuclear line emission. lEdwards et al.l 
( l2007h find that only BCGs at the centers of cooling core clusters 
are more likely to display emission lines than other massive (clus- 
ter) galaxies. They also confirm that in optically selected cluster 
samples (such as our sample), the BCGs are not more likely to dis- 
play emission lines than massive comparison galaxies. 

Some of the line emission could be attributed to cluster cool- 



ing fl ows instead of optical AGN activity (e.g. IVoit & Donahue! 
1 19971) . In addition, only half of the galaxies have detectable 
emission lines, so we would prefer a stellar age indicator that 
can be measured for all galaxies. The strength of the Balmer 
break £'„(4000) is measurable with a high SNR in all the galaxy 
spec tra and is an indication of the age of the stellar popula- 
tion. iKauffmann et al ] l l2003al) find that galaxies separate into two 
distinct populations, with young, star-forming galaxies having 
£'„(4000) < 1.6 . We find that the number of galaxies with 
£)„(4000) < 1.6 is very similar to the number of star-forming 
galaxies identified from the BPT diagram, for both the BCGs and 
non-BCGs (see Fig. In Fig. [TH we plot Z>„(4000) against 
Mu - Mg within the zso23 aperture. M„ - Mg also straddles the 
Balmer break and can thus serve to probe the average stellar pop- 
ulation at radii larger than the fiber aperture. The percentage of 
galaxies with blue M„ - Mg overall colors is compatible for BCGs 
and non-BCGs, i.e. again there is no indication for enhanced star 
formation in BCGs. Since we have matched in Mg - M,. color and 
have argued that this is essentially a match in stellar mass-to-light 
ratio, we do not expect a systematic difference in M„ - Mg (and thus 
this exercise may serve to confirm this argumentation). Without the 
match in Mg - Mr, the comparison sample contains more spiral 
galaxies, i.e. more star-forming galaxies than the BCG sample. 

To conclude, we do not find evidence for increased amounts of 
star formation in BCGs. However, we would like to caution that our 
sample is not well suited for such an investigation: we probe only 
the very center of the galaxies (and thus cannot detect line emission 
on larger scales), and we do not have X-ray data for our sample 
(previous studies suggest there is a strong connection between the 
cluster X-ray properties and star formation). 
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7 SUMMARY AND DISCUSSION 

We have developed a refined algorithm for selecting brightest clus- 
ter galaxies (BCGs) in the C4 cluster catalog, and we have im- 
proved the determination of velocity dispersion and cluster mem- 
bership. This refined cluster sample consists of 625 galaxy clus- 
ters at z < 0.1, and spans a wide range in velocity dispersion, 
from galaxy groups to rich clusters of galaxies. This, along with 
the detailed information available from SDSS for member galaxies, 
makes it a suitable local comparison sample for optically-selected, 
high-redshift cluster samples. 

Since the original SDSS magnitude measurements of BCGs 
are affected by excessive sky subtraction, we have developed a pro- 
cedure to recover more accurate magnitudes by adding a fraction 
of the difference between the local and global sky background 
estimates to the radial surface brightness profiles of the SDSS 
galaxies, and by determining new magnitudes from these corrected 
profiles. We show that this procedure removes the systematic bias 
in z - y color as a function of angular size for elliptical galaxies 
with photometry from both the SDSS and the 2MASS surveys. 
We also show that our reconstructed surface brightness profiles of 
BCGs agree well with previously published aperture photometry, 
at least to the isophotal limit of r = 23mag/n" within which we 
measure the flux. 

The properties of BCGs are governed by two main factors: 
their large (stellar) masses and their locations at the bottom of the 
potential well of their host cluster. Our large sample of BCGs en- 
ables us to disentangle the influence of these two factors and to 
assess the extent to which BCGs differ from 'ordinary' massive 
galaxies. 

We investigate the occurence of radio-loud AGNs in the BCG 
population and we find that BCGs are more likely to be radio- 
loud than other galaxies of the same stellar mass. This enhance- 
ment ranges from a factor 10 at stellar masses of 5 x IO'^Mq to 
less than a factor of two at masses larger than 4 x 10" Mq. This 
difference is arguably the most striking difference between BCGs 
and non-BCGs, and is likely due to the additional hot gas available 
in groups and clusters to fuel the radio AGN. 

The influence of the cluster environment is also evident as an 
increase in the fraction of dark matter in BCGs. The main obser- 
vational signature of this excess dark matter is that BCGs have 
larger radii and higher velocity dispersions than non-BCGs. How- 
ever, since the conversion of these quantities to dynamical mass 
depends on the shape of the mass and light profiles of the galaxy, 
this result could also be mimicked by non-homology between the 
BCGs and non-BCGs. Either case leads to a different slope of the 
fundamental plane for BCGs, one that is much closer to the virial 
plane than the observed fundamental plane of normal ellipticals. 

It is interesting to note that the differences from 'normal' 
galaxies are particularly evident in lower-mass BCGs: these are the 
galaxies that have a factor 10 higher probability of being radio-loud 
and have a significantly larger dark matter mass to light ratio (or a 
larger degree of non-homology) when compared to the non-BCGs. 
The low mass BCGs also deviate most strongly from the generic 
fundamental plane. 

We find that the slope of the Faber-Jackson relation is differ- 
ent for BCGs, in that their velocity dispersion rises less steeply for 
a given increase in luminosity than for non-BCGs. We also find evi- 
dence that this effect is stronger for BCGs in massive clusters. Such 
a change in the Faber-Jackson relation is predicted if these systems 
form in dissipationless mergers along elliptical orbits. Our results 



thus support the scenario where BCGs form mainly via dissipation- 
less mergers, and imply that the merger orbits are preferentially ra- 
dial in the most massive clusters. 

A difference in the Faber-Jackson relation also implies that 
BCGs can follow at most one of the power-law relations often 
used to estimate the mass of the supermassive black hole at their 
center, i.e. either the Mbh - cTv or the Mbh - L relation. In fact, 
this is already obvious from the parameter distributions shown 
in Figs. [To] and [TT] for the same distributions in stellar mass and 
in luminosity, we find a systematically different distribution in 
velocity dispersion for our BCGs. The standard re lations between 
black hole mass and bulge velocity dispersion jTremaine et al.l 
120021) have likely been derived for non-BCGs and may lead to 
systematically wrong black hole estimates for the brightest group 
and cluster galaxies. 

We find that BCGs have very similar mean stellar ages and 
metallicities to non-BCGs. They have slightly higher a/Fe ratios, 
indicating that their stars may have formed over a shorter time inter- 
val. Finally, BCGs display weaker optical emission lines than non- 
BCGs of the same stellar mass. In both BCGs and non-BCGs, the 
detected emission lines stem predominantly from low-lumin osity 
optical AGNs. In the accompanying paper jBest et alj|200^ , we 
further investigate the occurence of AGN activity in BCGs, and 
argue that the radio-loud and the emission-line AGN activity are 
independent, unrelated phenomena. 
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We thus need a 'trigger' to determine for which objects the 
local sky should be kept, and for which it should be replaced by 
the global sky. 

We make the trigger for each galaxy / a function of the ra- 
tio LR, of its luminosity L, to the luminosity L„b of its neighbors 
(within r.6), defined in the following way; 



LR; 



-2.5 log(L, /L„b) 

j e {sjalaxiesi jE (stars) 

— + 10 



j e IgidaxiesI 



2 

j e {stars) 



APPENDIX A: IMPROVING SDSS MAGNITUDES 



As described in Sect. 12.21 we need to correct the photometry of the 
BCGs and around 200,000 galaxies at z < 0. 1 for the fact that the 
standard SDSS photometric pipeline overestimates the sky back- 
ground. 

The SDSS photometric pipeline photo estimates the global 
sky within a frame (2048 x 1498 pixels; 13'.5 x 9'.8) from the me- 
dian value of the pixels in that frame, clipped at 2.3cr. The local 
sky background is then determined with the same sigma-clipping 
within a box of 256 x 256 pixels (r.7 x 1'.7) on a grid every 128 
pixels, and interpolated between these positions. This sky estimate 
is then subtracted from the image, and the photometry is performed 
on the sky-subtracted image. 

If a large fraction of the pixels in a 256 x 256 pixels box is part 
of an object rather than blank sky, this procedure causes the local 
sky to be an overestimate of the true sky background. This may 
happen for large galaxies, in crowded fields, and also around stars, 
since the wings of the PSF result in a considerable stellar halo. 

An overestimation of the sky background results in an under- 
estimate of the surface brightness of the object's pixels - thus, the 
total effect this has on the flux of an object scales with the square of 
its radius. T his effect is therefor e particularly severe for large galax- 
ies. Indeed, Lauer et al] 1 2006 ) anal yze the photometry of BCGs 
analyzed in jPostman & LaueJ ( ) 19*95!) that lie in the SDSS DR4 and 
show that the discrepancy in the estimated luminosity is a function 
of BCG radius (see also below). 

It has been suggested that the true flux of an object may be re- 
covered by using the global sky estimate instead of the local 
sky estimate (Masataka Fukugita; SDSS mailing list). This ap- 
proach should allow a new flux measurement without performing a 
new, independent photometric analysis of the raw images. For each 
object, the ID radial surface brightness profile, measured in 15 ra- 
dial bins, is available. By adding the difference between local sky 
and global sky to these surface brightnesses, the flux of an object, 
assuming the global sky level in the respective field, can be mea- 
sured. 



where r,, , is the petrosian radius of the galaxy and dj j is the dis- 
tance between galaxy ; and its neighbor j. Thus, the contribution of 
a neighbor to L„b is weighted by a Gaussian of width equal to twice 
the petrosian radius of the galaxy / - this is the aperture within 
which the SDSS petrosian flux is measured. Since we find that the 
presence of a star close to the object results in a substantial over- 
estimate of the galaxy luminosity if the global sky is used, stars 
are weighted with an additional factor of 10 (the exact value of this 
factor makes only little difference, but we do find slightly better 
results using a factor of 10 rather than 5). We suspect that the stel- 
lar halo due to the broad PSF wings is accounted for primarily as 
local sky background. 

We find that for LR, < -2.5, the flux of neighboring objects 
is negligible enough for the global sky to be the superior sky 
estimate. For LR, > -1, the flux of neighboring objects contributes 
a substantial fraction of the local sky estimate, so that it cannot 
be substituted by the global sky. 



A2 The method 

Rather than simply substituting the global sky for the local sky 
(or not), we calculate the fraction fsky of the difference between 
local and global sky to be added to the surface brightness profile 
of each galaxy according to the following criteria: 

(i) If Asky = skyioj.j,i - skygiobai < in any one of the five bands, 
then /,ky = . In these cases, the two sky measurements are 
essentially equivalent, and subtracting flux from the surface 
brightness profile will add noise to the flux measurement. 

(ii) Values of As^y are constrained to be < 10"'maggies/n". 

(iii) The maximum value for /^ky is 0.7 . In Sect. lA4l we show that 
this is superior to using /JJ^" = 1 . This value is assigned to 
objects with LR, < -2.5. 

(iv) The minimum (non-zero) value for /™ is 0. 1 . This value is 
assigned to objects with LR, > - I. 



(v) For objects with -2.5 < LR, < -1, /si^y is a linear function of 



LR,, being continuous at the endpoints with /J^" and /J 



rmin 
^sky 



Al Neighboring objects 

For isolated objects, the above argument implies that the global 
sky is generally a better sky estimate than the local one. However, 
for blended objects, we find that the local sky accounts for a large 
fraction of the flux from the respective neighbors; in these cases, 
therefore, using the global sky estimates would include flux from 
the neighbors and thus lead to an overestimate of the luminosity 
of the object. This is particularly true if the neighbor is a star of 
similar brightness to the object in question. 



A3 Comparison for the BCGs of lPostman & Lauerl jl995l) 

In order to assess the performance of our method to correct the 
SDSS magnitudes, we need (an) external dataset(s) with accurate 
photometry. We rely on t wo such datasets: t h e aper ture photometry 
of BCGs published in IPostman & Lauei) (ll995L referred to a s 
PL95 hereafter), and the 2MASS survey jSkrutskie et"ai]|2006h . 
The method is considered successful if it can simultaneously 
reproduce the curve-of-growth of a large fraction of the PL95 
BCGs, and if galaxies do not show any systematic bias in their 
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median SDSS-2MASS colors. 

In Fig. lAll we demonstrate that our method can reproduce the 
aperture photometry of the PL95 BCG sample. Of the 1 19 BCGs in 
that sample, 35 have imaging data from DR5, and 12 are also in our 
BCG sample (which was based on DR3). We show the curve-of- 
growth of six representative BCGs in Fig. lAll We as sume r' -Rg = 
0.25 as the typical color of an elliptical at z = jpukugita et al.l 
[1991 in order to compare measurements taken in different bands. 

The examples show how the original SDSS photometry breaks 
down at radii larger than ~ 20". Our improved photometry is able 
to reproduce the PL95 aperture photometry much better, out to 
> 80". Of the 12 BCGs in common with our sample, our pho- 
tometry fails to accurately reproduce PL95 only for the BCG of 
Abell 160 (C4_DR3.2025): this is because PHOTO attributes a large 
fraction of the BCG luminosity to a secondary nucleus; hence LR 
is rather large, and only little sky background is added to the bright- 
ness profile. Of the 35 BCGs in DR5, this is also the case for the 
BCG of Abell 1 185 (not shown). The crowding in the fields around 
the BCGs of Abell 2040 and Abell 2052 also leads to an underesti- 
mation of the flux by about 0.2 mag (not shown). 

In addition to bad sky subtraction, noise in the surface bright- 
ness profile plays an important role particularly at large radii, as can 
be seen in the curve-of-growth of the BCG of Abell 779 (Fig. lAU . 
The last non-zero bin of the SDSS surface brightness profile causes 
a jump in total magnitude of more than 0.5 mag, which is not seen 
in the PL95 photometry. To avoid such jumps, we compare the cu- 
mulative flux for each bin of radius > 18" to the flux predicted 
for that bin from the slope of the curve-of-growth from the two 
bins prior in radius. If the flux is more than 30% larger than pre- 
dicted (or if the surface brightness in the respective bin is negative), 
the measured flux is substituted by the predicted flux. For several 
BCGs this procedure improves the agreement between the SDSS 
and PL95 brightness profiles. For one (the BCG of Abell 2147), 
this causes disagreement at large radii. However, this occurs at a 
significantly larger radius than the one within which we measure 
the magnitudes. 



A4 Comparison to 2MASS magnitudes 

Of the 200,000 unique SDSS galaxies at z < 0.1, about half have 
a counterpart in the 2MASS survey's Extended Source Catalog 
(XSC); we made use of the list of 2MA SS objects matched to SDSS 
galaxies provided via the NYU VAGC l lBlanton et alj2005l) to iden- 
tify these galaxies and obtain their 2MASS properties. 

In order to be able to directly compare the SDSS and XSC 
data, we work with isophotal magnitudes. One of the magni- 
tude measurements provided by 2MASS is measured within the 
K = 20mag/n" isophote. A" = 20 corresponds approximately 
to r = 23, hence we measure SDSS magnitudes within the ra- 
dius where the surface brightness profile drops to r = 23 mag/n". 
We compare these magnitudes for elliptical galaxies, selected us- 
ing standard cuts that have been used in previous SDSS studies 
(i.e. CsDss = ^^o°''/'^?;)"' > 2-86 and g^^l^ - r^DSS). We also limit 
the sample to galaxies where the SDSS and 2MASS radii agree to 
within 30% (this criterion changes for each magnitude measure- 
ment). We consider only galaxies with a radius larger than 7", the 
minimum radius n ecessary to avoid PSF effects for 2MASS data 
jjarrett et al]|2000l) . After correcting for galactic extinction (and 
converting the 7-band magnitude to an AB magnitude), we com- 
pare the SDSS z-band magnitudes to the 2MASS /-band magni- 
tudes, since these bands are adjacent in wavelength and thus color 




isophotal radius ["] 

Figure A2. The difference between SDSS z. magnitude (measured within 
the r = 23 mag/n" radius) and 2MASS J magnitude (measured within the 
K = 20 mag/n" radius) as a function of galaxy size (given by the r = 
23mag/n" isophotal radius). We measure the median difference in bins of 
radius; the error bars denote the 68% levels of the distribution in each bin, 
divided by the square root of the number of galaxies in that bin. For the 
black symbols, the sky subtraction has not been changed from the original 
SDSS value. For the red symbols. Z^™" = 1-0, i.e. if an object meets the 
criteria to replace the local by the global sky. the full difference is added 
(see Section lA2l for details). For the green symbols, f"™ = 0.5, and for the 



blue symbols, /J™" : 
values for all galaxies. 



0.7 and/™" : 



0.1. The dashed lines show the median 



differences from different stellar populations can be expected to be 
minimal. 

In Fig. lA2l we investigate the color (z,23 - /a:2o) as a function 
of galaxy size (given by the r = 23 mag/n" isophotal radius in the 
SDSS) for four different sky subtractions. Clearly, for the original 
SDSS sky subtraction, the color term is a strong function of galaxy 
size, indicating that the sky is systematically oversubtracted in the 
optical band. This accounts for a systematic difference of the order 
of 0.1 mag for galaxies larger than 20" (for i ndividual galaxies , 
this may be much more, as demonstrated by iLauer et al]|2006t) . 
However, substituting the local sky estimate by the global one 
(i.e. /J^y" = 1.0) leads to an underestimation of the sky, and thus 
an overestimation of the luminosity. We find that with /J^" = 0.7, 
there is little variation of the median color term with galaxy size, 
and no systematic trend (note that setting = 0. 1 makes only a 
small difference). 



A5 Final magnitudes 

The isophotal magnitudes are not redshift independent because of 
cosmological surface brightness dimming. We thus modify the al- 
gorithm to correct both for (1 -I- z)* surface brightness dimming and 
for galactic extinction (i.e. these corrections applied directly to the 
radial profile before the magnitudes are measured). 

We choose to use isophotal magnitudes corresponding to r = 
23 at z = (i.e. r = 23.41 at z = 0. 1). This is a rather bright cut-off. 
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Abell 1 139, C4 1 142 
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Figure Al. Comparison of the aperture photometry of jPostman & Laue (shown as diamonds) and SDSS photometry for 6 BCGs. The uncon'ected 

SDSS magnitudes are shown as triangles, with the respective cubic spline fit as a dashed line. By adding flux of /skyAji^y, we obtain the curve-of-growth shown 
by + symbols and as dash-dotted line. Noisy photometry at large radii can cause jumps in the curve-of-growth; this can be alleviated by the method described 
in the text, and is shown as X symbols, and as solid line. The dotted lines indicate the radii at which the surface brightness profile reaches a magnitude of 
r = 23 mag/n" (left line) and r = 24 mag/n" (right line). 



but is less sensitive to residual errors in the sky subtraction and/or 
surface brightness measurements (cf. Fig. lAlt . 



A6 Influence of the sky subtraction on our results 

Our analysis of BCG properties (Sect. [6]( as a function of stellar 
mass depends critically on the correct determination of the lumi- 
nosities of the objects in our sample. If the luminosities we assign to 
BCGs are underestimated, then the comparison galaxies (which are 
matched in stellar mass) would be systematically less massive than 
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Table Al. Comparison of the results of Kolmogorov-Smirnov tests (expressed as the decimal logarithm of 1 minus the confidence level at which the null 
hypothesis, that the distributions are drawn from an identical parent population, can be rejected) on BCG and non-BCG distributions of various parameters, 
for different versions of sky subtraction applied to the SDSS photometry. 



sky subtraction 
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-14.46 
-12.76 
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-2.61 
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-2.27 
-2.19 


-14.00 
-12.57 


-0.57 
-0.61 


-2.71 
-2.74 


-1.04 
-1.75 


-0.04 
-0.01 



the BCG, leading to possibly spurious differences in their physical 
properties. 

To estimate the effect that sky subtraction has on our results, 
we repeat the analysis presented in Sect.[6]with photometry derived 
with values of = 0.5 and /™ = 1.0. To evaluate whether there 
could be significant quantitative effects on our results, we compare 
the results of the Kolmogorov-Smirnov test on the BCG and non- 
BCG distributions for several parameters (as shown also in Fig- 
ures [To] and lilt in Table lAll We find only minor differences to 
our initial results; none of our conclusions are altered as a result of 
using a different value of X,™"". 



APPENDIX B: CLUSTER EXAMPLES 

In the following we present finding charts and redshift histograms 
for clusters referred to in the text. The finding charts are centered 
on the BCG, and spectroscopically confirmed cluster members are 
marked by triangles. In the redshift histograms, a short-dashed line 
indicates the cluster redshift, the long-dashed line the redshift of the 
BCG, and the dotted lines the ±3cr,, ^.i limits. Also, the redshift (both 
the C4 and our measurement) as well as the velocity dispersion 
measurements (given in km/s) are listed in the redshift histograms. 
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z = 0.058 
z„ = 0.058 
a. = 627*11 
(J..C, = 738 



Figure Bl. C4_2003. Since most galaxies seem to cluster around the 
galaxy shown at the center of this image, this is the galaxy we identify 
as the BCG, even though the galaxy about 6' west-south- west of the 
center (marked by two lines) is brighter by a third of a magnitude. The 
R200 of this cluster is 1.5 Mpc, which translates to ~ 20'. 
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C4 1283 
z = 0.095 
z„ = 0.098 
a, = 754*55 
(J.,,, = 909 



101 



Figure B2. C4JDR3.1283. For this cluster, tlie C4 mean galaxy 
(marked by diagonal lines in the finding chart) lies 2.8 Mpc from the 
BCG which we identify. The large-scale distribution of galaxies at the 
cluster redshift suggests that C4 selected part of the infall region of 
this cluster, but fails to pick up the cluster itself. Note that the brightest 
galaxy in the field is IC 0504 at ; = 0.013, whereas many of the other 
bright galaxies in the foreground belong to C4_DR3_1356 at z = 0.03. 
Below, we show thumbnail images of the C4 mean galaxy (left) and our 
BCG (right). 
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C4 2163 
z = 0.070 
z„ = 0.079 
t'. = 225-^, 
= 2289 



Figure B3. C4.DR3 J163. In the C4 catalog, this cluster is Usted as 
a structure at z = 0.079, with o"v_ci > 2000 km/s. The mean galaxy is 
marked in the finding chart below and is part of the group of galaxies 
at z = 0.07, which our algorithm identifies. In the redshift histogram, 
we show the distribution of galaxies within projections of IR200 (black, 
shaded), 3i?2oo (red), and 5i?200 (blue). Obviously, the difference in the 
redshift measurements as well as the large C4 velocity dispersion are 
due to another structure at z ~ 0.082. In Fig. IB4I we show that this 
background structure is associated with C4_2124. 
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Figure B4. C4_2124. The finding cliart below is not centered on tiie 
BCG of C4.2124, but on the BCG of C4J3R3.2163, in order to illus- 
trate the sheet-like structure at z ~ 0.082 which led to the deviating red- 
shift measurement for C4.DR3.2163. The BCG of C4.2124 is marked 
in the finding chart, and its S200 (equivalent to 1 Mpc) is indicated by a 
dashed circle. 
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